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G^ ! Abstract 

Models generalizing the su{2) XX spin-chain were recently introduced. These XXC mod- 
O I els also have an underlying su{2) structure. Their construction method is shown to generalize 

■ to the chains based on the fundamental representations of the Am Lie algebras. Integrability 

ly^ I of the new models is shown in the context of the quantum inverse scattering method. Their 

i?-matrix is found and shown to yield a representation of the Hecke algebra. The diagonal- 
ization of the transfer matrices is carried out using the algebraic Bethe Ansatz. I comment 
. on eventual generalizations and possible links to reaction-diffusion processes. 
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1 Introduction 



In the course of generalizing the Hubbard model new models were discovered: the so-called 
XXC models ^, These one-dimensional (spin-chain) integrable models have a natural 

expression in terms of su{n) generators, rather than higher-spin representations of su{2). The 
first hamiltonian of the integrable hierarchy is bilinear in terms of su{n) generators. The sym- 
metries of the models are generators of su(ni) su{n2) © ^^(1)- The algebraic Bethe Ansatz 
diagonalization of the transfer matrices requires a nesting similar to the XXZ models based on 
the fundamental representation of su{p). Finally, for a subclass of the XXC models the explicit 
form of all the conserved charges was found [Q]; their expressions in terms of su{n) structure con- 
stants indicate an su{n) interpretation. New identities for the structure constants were derived 
as a by-product. These should admit further generalizations. 

However the XXC models also share features associated with the spin-^ XXZ model and 
appear as a kind of higher-dimensional representations of the i?-matrix of the spin-^ model. 
Indeed, the latter model is a special case of the XXC models and the /^-matrices of these models 
share a common structure, with their building blocks satisfying the same algebraic relations. 
Moreover the algebraic Bethe Ansatz can also be interpreted as being nest-less and therefore a 
simple generalization of the su{2) one. Another argument in favor of the su{2) interpretation 
can be found in [||. One obtains the quadratic hamiltonian of a particular XXC model, with 
open boundary conditions, as the 'infinite-coupling' restriction of the Hubbard hamiltonian on 
a subspace of the complete Hilbert space. The authors of ||5| showed that this model possessed 
a surprisingly large affine symmetry based on su{2). A generalization of the resulting model 
yielded a subclass of open-boundaries XXC hamiltonians; however their symmetries favor an 
su{n) interpretation. 

One therefore may try to generalize the construction method used in deriving the XXC 
models, to the XXZ models which are based on the fundamental representations of the Am Lie 
algebras. This approach turns out to work and the resulting models are obtained and studied 
in this paper. One starts with the Am i?-matrices. Their structure allows a straightforward 
generalization retaining their operatorial form and their Am characteristics. The new matrices 
are given and shown to satisfy the Yang-Baxter equation. Integrability of the models is then a 
simple consequence of the quantum inverse scattering framework. The symmetries are obtained 
and the transfer matrices are diagonalized by algebraic Bethe Ansatz. I conclude with some 
remarks and possible physical applications to reaction-diffusion processes. 



2 New models 

The Yang-Baxter equation (YBE) is at the center of Quantum Inverse Scattering Method 
(QISM) used to obtain quantum integrable one-dimensional spin-chains and their covering two- 
dimensional classical statistical models ^, 0, |9|. There are now several methods which can be 
used to obtain solutions, R- or L-matrices, of the YBE. An important and quite general method 
relies on the use of affine quantum groups based on Lie algebras. Rather than directly solve 
the cubic equations resulting from the YBE, one solves linear equations where R appears as 
the intertwiner between two possible deformed coproducts (tensor products). This method was 



used in particular in [10| and explicit trigonometric solutions were obtained for the fundamental 
representations of the classical Lie algebras Am, Bm, Cm, Dm, and their twisted versions. 

For the untwisted algebra Am-i, the trigonometric ^-matrix of the fundamental represen- 
tation is m^-dimensional and can be found in |To|] : 

R{y) = sin(7)(y^£;^^®E"" + 2/-i^^'3^®E"") 

a>f3 a<l3 
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+ sin(A + 7) + sm(A) ^ O E"'^ (1) 

where y = e*^, A is the spectral parameter and 7 the quantum deformation parameter. The 
E'^^ are m x m matrices with a one at row a and column /? and zeros otherwise. The matrix R 
satisfies the Yang-Baxter equation 

Al2(A)^23(A + ^i)RM = i?23(/x)^12(A + /i)i?23(A) (2) 

for any fixed value of 7. Here and in (|l2|) , the notation Oij {i ^ j) means that the operator 
O acts non-trivially on the i^^ and j^^ spaces, and as the identity on the other spaces. The 
regularity and unitarity properties also hold: 

i?(0)=Isin7, i2(A)i?(-A) = I sin(7 + A) sin(7 - A) (3) 

I now give new solutions to the YBE which are obtained by a multi-state generalization of 
expression (|l|). One first rewrites the latter matrix as 

R(X) = (yp(+) + sin7 + sin(A + 7) + P^^^ sin A (4) 

and looks for representations of the operators P which allow R to satisfy the Yang-Baxter 
equation. Let rij be m positive integers such that 

m 

rii = n and 1 < ?^l < ... < ?^m < n — 1 (5) 

i=l 

The inequality restrictions avoid multiple counting of models. Split the set of n basis states into 
m disjoint sets Ai: 

card (Ai) = Ui , A n ^ j = for i / j (6) 
Ai should not be confused with the Lie algebra su{i + 1). Consider the following expression for 

p(3). 

The twist parameters Xa^a- are arbitrary complex numbers. The remaining operators are given 
by: 

p(i) _ (p(3))2 ^ p(+) ^ p(-) = H (-^"'°' ® + ® ■^"'°') (8) 

l<i<j<mai£Ai O-j&Aj 

m 

p(2) = I _ pw = Y, Y ® (9) 

The operators and P^^^ correspond respectively to the sums over the first and second 

summands in (^) . There is also another way of writing (Q) : 

R{X) = I sin(A + 7) +P sinA (10) 
P = P(3) _ (g-i7p(+) + gnp(-)) (11) 

A straightforward if tedious calculation shows that: 

P2 = -2 P cos 7 , P12P23P12 + P23 = P23P12P23 + P12 (12) 
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P is therefore a generator of the Hecke algebra. These relations imply that the Yang-Baxter 
equation is satisfied. The regularity and unitarity properties @j still hold. I denote this model 
by {ni,...,nm;m,n). 

There is a simple graphic and mnemonic representation of the foregoing operators. To each 
state assign a point in the plane. States belonging to the same set Ai are not linked while those 
belonging to different sets are linked. A given link corresponds to a given summand appearing 
in the expression of P^^\ a tensor product of two step operators, and also to the summands 
appearing in P^^\ a tensor product of diagonal operators. Similarly the summand of P^"^^ 
corresponds to missing links in the diagram. Links and missing links exhaust all possible links 
which could be drawn between the n states. This representation is illustrated with two examples 
in figure 1. One can also read this diagram as follows. One starts with an su{m) system and 
replaces every state with an arbitrary number of copies. The copies originating from the same 
state do not 'interact' among each other; they interact with all other states and their copies as 
dictated by the original diagram. 




Fig.l: The diagram on the left corresponds to the system (2, 3; 2, 5) and the one on the right to 
(1, 2, 2, 3; 4, 8). Here Ai stands for A- 

The choice of which states go into which set Ai does not yield inequivalent systems. For 
instance the two systems, {Ai = {1} , A2 = {2,3}) and {Ai = {2}, A2 = {1,3}), are related 
by a simple (orthogonal) permutation matrix whose A^-fold tensor product with itself yields the 
unitary matrix which relates the two A^-sites integrable models 

The operators P^^^ and P^^^ form a complete set of projectors on the tensor product space 
C" ® C": 

p(l) _^ p(2) ^ J ^ (^p(l))2 ^ p(l) ^ ^p(2)^)2 ^ p{2) ^ p(l)p(2) ^ p(2)p(l) ^ q ^-^^-^ 

One also has (p(^))^ = P^^\ However, for m > 2, the operator P^'^' does not satisfy the 3-sites 
relations of the 'free-fermions' algebra A found in Q. This is an important difference between 
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m = 2 and m > 2. In the latter case no 'conjugation matrix' exists and it does not seem possible 
to couple two such models. 

For m = 2 one obtains the XXC models in their asymmetric guise, as there are factors of 
y and y~^. The transformation to the symmetric models of is given by a simple 'gauge' 
transformation: 

rGT^X) = {A{X) 1) R{X) (I ^(-A)) (14) 
where A{\) = J2^_^ ^"laig^^ci _^ ^a2a2giAc2 ^i^-j^ C2 - ci = 1. All the above properties are 
preserved by such a transformation. For m > 2 however it is not possible to remove the y^^ 
factors. For m = n and all parameters x equal to one, one obtains the matrix (|l]). 

The rational limit of R is obtained by letting A 7A, dividing by sin 7 and taking the limit 
7 — 0. These manipulations conserve all the properties of the i?-matrix. In particular, one 
obtains 

R{X) = p(i) + (1 + A)P(2) + AP(3) (15) 
i2(0) = I , R{\)R{-\) = I (1 - A2) (16) 

If all parameters Xa^aj are equal to each other, both the trigonometric and rational matrix 
have the following symmetry 

m m 

[M ® M, R{X)] = , M = J2 M^^^ = H "^Sk^"'^' 

k=l k=l afc,/3fce.4fc 

If however two blocks M^'^'^ and M^'^ \ k < k , are diagonal M is a symmetry of R with the 
corresponding parameters Xa^a , being unconstrained. 

The transfer matrix is the generating functional of the infinite set of conserved quantities. 
Its construction in the framework of the Quantum Inverse Scattering Method is well known. 
The Lax operator on a chain at site i with inhomogeneity /Xj is given by: 

Loi{X) = Roi{X - fii) = Voi ^oi(A - fii) (18) 

where V is the permutation operator on C" iX) C". The monodromy matrix is a product of Lax 
operators 

r(A) = MoLojv(A)...Loi(A) (19) 

where is the number of sites on the chain and is the auxiliary space. The transfer matrix 
is the trace of the monodromy matrix over the auxiliary space: t(A) = Tro [T(A)]. The intro- 
duction of M corresponds to integrable periodic Af -twisted boundary conditions. A set of local 
conserved quantities is given by 

The YBE implies the following intertwining relations for the elements of the monodromy matrix: 

R{Xi - A2) r(Ai) r(A2) = r(A2) r(Ai) r{Xi - A2) (21) 

Taking the trace over the auxiliary spaces, and using the cyclicity property of the trace, one 
obtains [t(Ai), r(A2)] = 0. The hamiltonians Hp therefore mutually commute. 

The quadratic hamiltonian calculated from ([20|), for /Xj = and M = I, is equal to 

H2 = Y: Hn^i = E (^^1 + P^ii 7 + (^;i -P^Jt^)^ sm 7) (22) 



4 



For \xa/3\ = 1 and 7 purely imaginary the hamiltonian sin 7 x H2 is hermitian. For m = 2 this 
hamiltonian is also hermitian for real values of 7. Under the periodic boundary conditions the 
non-hermitian part does not contribute. This is easily seen from the transformation ( p!4[ ) and the 

fact that the hamiltonian density -ffjj+i is equal to the derivative at zero of R{\). The rational 

(o\ ('2') 

limit yields: H2 = + Pjj+i)^ and, provided \xa(}\ = 1, the hamiltonians are hermitian. 

The cubic conserved quantity is obtained from ([20|) by a direct calculation. One finds = 



— i-^j-iji -^jj+i] ~ ih?^^' '^here Hjj^i is the hamiltonian density of (22) or its rational 
version. The sin^ 7 is replaced by 1 for the rational limit. The commutator can be easily derived 
usmff 

= Sp^E°'^. No general closed form expressions for the higher conserved quantities 
have yet been derived in the literature; however see for some specific cases, and the references 
therein for related issues. 

The transfer matrix, and therefore all the conserved quantities, have the symmetries of the 
^-matrix. Define the magnetic- field operators as: 

jjaA^j2j^aA ^ ak,PkeAk (23) 

i 

One has the following commutation relations for both the trigonometric and rational forms: 

, r(A)] =0 if and only if (24) 
yj < k, \/-fj € Aj x-y^ak = x-y^f3k and Mj > k, M-fj E Aj x^^^^ = xp^^^ 
and nialak = mf^f^^ , m^^k =0 ^Ik 7^ at , mf^^^ =0 \f-/k 7^ 

In particular the diagonal operators commute with the transfer matrix without any 

constraint on the twist parameters, but with the above constraints on the matrix M. The 
rational transfer matrix, with all Xq-q,^. 's equal to one, may have additional symmetries. For 
k ^ k' one finds 

^jj'^^^u' ^ ^rat.(^)] ^ [i^f*'"", r>^^t-(A)] = (25) 
if and only if = n^/ = 1 and m'alak = ^13 , 

k k 

In the following we shall concentrate on the case where M = I and all the x's are equal to 
a single parameter x, despite the fact some results hold for generic parameters. The full local 
symmetry then is su(ni) © ... © su{nm) © uil) © ... © u{l) where there are m — 1 n(l)'s. 



3 Algebraic Bethe Ansatz 

The diagonalization by algebraic Bethe Ansatz of the foregoing models combines features from 
the diagonalization of the (1, l;m,m) models and the XXC models. We refer the reader to 
1^, ^ |l^ for some details and give here the new features and results. The ^-matrix here is 
redefined as the matrix (^) divided by sin(A + 7). This adds —A^ cot 7 1 to the hamiltonian (|2^ ) 
and + (-0^ to H^. For the rational limit one adds — A^I and +NI respectively. 

Let k^^ € {!,... ,rn-} and 7 (i) S ^,(1) be given. The action of all the elements of the 

monodromy matrix on the pseudo- vacuum ||7,(i)) = |7,{i)) ■■■ ® |7,(i)) is easily derived: 

T^ W "fe 1 17,(1) ) 7^ , T „ I |7,{i) ) / (26) 
ftp 

Il7.«) , T,,..\W^^) = Ul,{x-^''^'-'o'^^^ (27) 



^7 (1)7 (1) IITJi)/ 
fe„ fc„ 
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\/k ^ and Va (i) 7^ 7 (i) . All other elements of T annihilate this vector. 



Let C/3j. = ,,^/3^.. Operator C/j^, exactly flips a state |7,(i)) into a state Thus these 

operators acting on ||7,(i)) give a linear combination of states where exactly one state in ||7,{i)) 
has been changed to \(3k), at every site. To show this one uses the following relations: 

[i?r"SC^/3^,] = Sa,f3^>Cp,, ^ VA:, VA:VC^ (28) 
^jj-yX' , C^J = 0, yk^ ki'^ , Va^(i) ^ 7,(1) (29) 



[^1 ° ^c-ft] = -Cft , yk^k',''. (30) 



'^0 

(1) 



Relation (pQ) also shows that (Ai)...C^(pi) (Ap^) ||7,(i)) = for > N. Relation ( ]29| ) shows 
that this same vector has no |a,(i)) state in it, and (|29| ) and (50) yield for a (i) 7^ 7,(i): 

Qi)(Ai)-C^(Pi)(ApJ IItjd) = 

The operators C^^, are therefore candidates for writing the eigenvector Ansatz. In contrast, the 

operator ^^^q, acting on ||7,(i)) changes the state of only the first site to |a^(i)). This is an 

kg kg 

unusual feature and these operators cannot be used to write down an eigenvector Ansatz. 
One may therefore take as Bethe Ansatz eigenvector 

|Ai,...,Ap,) ^F<''-<'c^(i)(Ai)...C^,p,)(ApJ ||7,(i)) (32) 

ki kp-^ 

where the parameters Aj and the coefficients F are to be determined. The sums run over all kt 
from 1 to m with ki 7^ /cg^'* and over 6^^ in Ak^- 

One then applies the transfer matrix on the state |Ai, Ap) and uses the algebraic relations 
(|2l|). The foregoing procedure, the nested algebraic Bethe Ansatz, is a cumbersome but straight- 
forward generalization of the one for the usual (1, 1; m, m) models, i.e. the su{m) XXZ model. 
The differences come from the sum on the multiple states in each set Ai, as already seen on 
the initial eigenvector Ansatz (|32[). The are m — 1 levels in the nesting and diagonalizing the 
transfer matrix at one level requires diagonalizing a new transfer matrix generated by the re- 
peated use of relations (21). This new transfer matrix corresponds to a system of the above 
type but with a reduced number of states and sites. The nesting stops at the last level at which 
the new transfer matrix is trivially diagonal. Technical considerations impose a decreasing or 
increasing sequence of one being needed for every level of the Ansatz. For the increasing 
sequence {k^^ = l,kl^^ = 2,...,m — 1), e = —1 in the eigenvalue and Bethe Ansatz equations 
while for the decreasing sequence e = +1. The sequence of systems appearing for the increasing 
and decreasing sequences are given by 

(ni, ...,nm;m, n) for pQ = N sites (n2, Um', m — l,n — ni) for pi sites — > 

• • • ^ {nm-i,nm; 2, n - ni nm-2) for Pm-2 sites (33) 

(ni, m, n) for pQ = N sites — > (ni, n^-i; w — 1, n — rim) for pi sites — > 

• • • ^ (ni, 77,2; 2, n - na rim) for Pm-2 sites (34) 
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The sequence of sites is non-increasing: 

iV = Po > Pi > • • • > Pm-l > (35) 

The Pm-i appears in the diagonahzation of the last system, (*, *; 2, *) in the above series; this 
system is an XXC one. The above procedure is akin to decimating the Am~i Dynkin diagram 
by going from either of its extremities to the other. 

The eigenvalue at one level is related to the eigenvalue at the following level as follows: 

)(A,{A«,...,A«}) = -i)n I " w ^ 



=1 Vsin(A-Af^+7), 
^'^/ x^sm{xt^''>^-X + ^) \ ^ ||^x-sin(A-Af))\ , 




t=i V sin(Af - A) ; ti Vsin(A - Af ^ + 7) / 

X A('"-'=-i'P*+i)(A,{Af+^\...,A(,^.+i)}) (36) 



sin(A - Af + 7) 



sin(A - Af +^)) 

(Ai°\ A^^^) = (^1, ...,/iAr) , k = 0,...,m-2 

where the subscript appearing in the (5-term is equal to A; + 1 (m — k) for e = — 1 (e = +1) 
respectively. The eigenvalue of the transfer matrix t(A) is A^"^'^'^\X, {Xj^\ \^^^}). The last 

eigenvalue A^^'P'^-^^X, {Xi" ^\ Xp^^_l^}) is independent the spectral parameter and of the 
inhomogeneities. It is an eigenvalue of the (constant) unit-shift operator sending the state on 
site i to site i + 1, on a lattice of Pm-i sites with rim (ni) states per site for e = — 1 (e = +1) 
respectively. 

The parameters A) appearing at every level are solutions of the Bethe Ansatz equations: 

sm(A,'«' - A,"') j,ii,.^si„(Af -A,'f-7)^ 




:t-<sm(Af' -A,'^.") 
=.n(Af' -At,"+7)^ 



and 




1, ...,pk , k = l, ...,m- 2 (37) 



. />(m— 1) ,(m— 1) , N 

sm(A,^ - + 7) 

(m— 1) ^(m— 1) 




(38) 

^{1) ^(pi) 

Finally, the coefficients F ''i ' ' are such that F and is an eigenvector of the transfer matrix 
t''"^~^'P^\X; Xi'\ Xp^). Note that, as usually happens in the ABA diagonahzation, these 
equations imply the vanishing of the residues of A('"~'''*''=)(A; {A^*^^}) at the Aj*^^^^'s. This was 
expected since the transfer matrix is non-singular at these values of the spectral parameter. 
All the possible combinations of piS satisfying ( |35| ) should be considered. For those with a 

first vanishing p^./ , equations (3^,37) truncate accordingly with A*-™'"'^ '^fc'-* = 1 and all products 
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n? set to one; equation ( pSD is not used. The BAE clearly display the Am-i Dynkin diagram 
structure when one uses the shifted parameters u^^^ = xf^'' + k'y/2. 

At every level / of this diagonalization one has the option of choosing among n (i) possible 

pseudo- vacuum. To cover the largest number of subspaces of the Hilbert space of the chain one 
should consider all choices. This gives distinct eigenvectors but equal eigenvalues. This reflects a 
large degeneracy of the spectrum and is different from the simple models (1, 1; m, m). Another 
difference lies in the appearance of the eigenvalue A^^'P^^^) at the last level. As explained in |^] 
the diagonalization of the unit-shift operator is in principle simple. Formulae for the degeneracies 
of its eigenvalues have been derived by M. Bauer |12|. 

Are there eigenstates not obtained by the foregoing procedure? From the action of the C 
operators one infers that states lying in the subspaces Ak iX) ... Ak, for every fixed k, are not 
'reached' by the Ansatz. I now fill a gap in and give the action of the transfer matrix on such 
states. One easily derives: 



r(A) |a«,...,arO 



N) 



(39) 



N 



+1 E n 



k'^k 



i=l 



sign(fc ~k) 



sin(A - fij) 
sin(A - /ij + 7) 



\a 



(1) 



k '-'"fc 



where q[.*^ G Ak and t^^*"'^^ is the unit-shift operator for states and N sites. Note that for a 
given k, all the states with the same eigenvalue for t^^'^'^) have the same eigenvalue for t(A). 

The rational limit of all the above equations is obtained by letting A* 7A*, fii — > 'yfii, 
dividing the eigenvalues by (sin 7)''^ and taking the limit 7 — > 0. 

For vanishing inhomogeneities, the N^^ power of the eigenvalue A('"'*''')(0, {0}) is equal to 
one because one has the unit-shift operator for sites and n states. One can verify this using 
equations (E 



The eigenvalues of the hamiltonians and higher conserved quantities are easily derived by tak- 
ing the logarithmic derivatives of the eigenvalue A('"'^'")(A, {a['^\ A^"*}) at vanishing spectral 
parameter. No closed form expressions of these derivatives are known. However, for vanish- 
ing inhomogeneities (/ij = 0), the first — 1 logarithmic derivatives of the eigenvalues in ( |3^ ) 
are easily found to vanish (in both trigonometric and rational cases). The first and second 
derivatives, E2 and E^, of the logarithm of the eigenvalue A^^^'P"^ can be easily found: 



and 



pi 



^2 = E 



sin7 



1 sin A^ ■* sin(Aj^'* + 7) 



pi 



1 



prat. \ ^ 

^' -feA«(A<ni) 



(40) 



PI ■ 

E - 2V sm7cosA,^^ 
' fesin2AS^)sin(Af)+7) 



pi 



i(A«)2(A« + l) 



{ET'f (41) 



4 Conclusion 

I have introduced new classes of solutions of the Yang-Baxter equation and found their symme- 
tries. These models appear as hybrids between su{n) XXZ and su{m) XXZ models but share 
their main characterizing features with the latter models. The diagonalization of the conserved 
quantities was then done using the algebraic Bethe Ansatz procedure. 

The issue of completeness of the Bethe Ansatz diagonalization, be it in its algebraic or 
coordinate form, is still an open issue for all but the simplest model, the 5^(2) XX model. 



8 



Various completeness derivations exists, however they all include some reasonable but unproven 
element. This should not be construed as a hindrance to the study of the thermodynamic limit 
in the Thermodynamic Bethe Ansatz (TBA) framework |13|. As the number of sites becomes 



large the Bethe Ansatz equations are transformed into a system of non- linear integral equations. 
To do this one assumes that the solutions of the BAE take a particular form for which their 
total number matches the dimensions of the Hilbert space. This form is true for most of the 
solutions and the results obtained are in complete agreement with other methods used to study 
the thermodynamic limit. This is probably due to the 'fact' that the set of irregular solutions 
has vanishing measure. 

Thus a detailed study of the spectrum in the TBA framework is desirable. Taking the 
logarithm of the Bethe Ansatz equations shows that the distribution of integers characterizing 
the solutions will get a contribution from the A(^'P'"-i) correction and degeneracies associated 
with the integers n^. Whether this influences the central charge and the conformal weights 
remains to be seen. 

Another open issue is the determination of a quantum group framework. A step in this 
direction was taken in |^ . It should admit generalizations and would shed some light on whether 
multistate generalization exist for higher representations of or for other Lie algebras. In 
particular, the multi-states Am models should extend straightforwardly to the Lie superalgebras 
su{m\n). 

Reaction-diffusion processes in one dimension are described by a time-dependent probability 
distribution for the configuration /?. This distribution obeys a stochastic master equa- 

tion which can be written as a Schrodinger equation with imaginary time, dt\P{t)) = —H\P{t)) 
and P{{P},t) are the components of the wave function in the basis of (species) states of the 
Hilbert space. The hamiltonian H contains the physical transition rates. For certain processes 
H was found to belong to integrable hierarchies described by a Hecke algebra [14, 15|. These 



hierarchies enter the class of models studied here. As the multi-states models are realizations 
of the Hecke algebra, possible physical applications may lie in the field of reaction-diffusion 
processes of multiple species. It would be interesting to pursue such an approach. 
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complicated eigenvalue degeneracy formula, and the Service de Physique Theorique of Saclay 
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